Leapover lengths and first passage time statistics for Levy fliglits 
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Exact results for the first passage time and leapover statistics of symmetric and one-sided Levy 
flights (LPs) are derived. LPs with stable index a are shown to have leapover lengths, that are 
asymptotically power-law distributed with index a for one-sided LPs and, surprisingly, with index 
a/2 for symmetric LPs. The first passage time distribution scales like a power-law with index 1/2 as 
required by the Sparre Andersen theorem for symmetric LPs, whereas one-sided LPs have a narrow 
distribution of first passage times. The exact analytic results are confirmed by extensive simulations. 
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The statistics of first passage times is a classical con- 
cept to quantify processes, in which it is of interest when 
the dynamic variable crosses a given threshold value for 
the first time, e.g., when a tracer in some aquifer reaches 
a certain probe position, two molecules meet to form 
a chemical bond, animals search for sparse food loca- 
tions, or a share at the stock market crosses a preset 
market value P, 0, Here, we revisit the first pas- 
sage time problem for processes with non-trivial jump 
length distributions, namely. Levy flights (LFs) and de- 
rive exact asymptotic expressions for the first passage 
time density Pfij) of symmetric and one-sided LFs. For 
the former, we obtain the Sparre Andersen universality 
Pfir) ~ T~^/^, while a narrow behavior is found for one- 
sided LFs. Apart from calculating the first passage times, 
we investigate the behavior of the first passage leapovers, 
that is, the distance the random walker overshoots the 
threshold value d in a single jump (see Fig. [1]). Surpris- 
ingly, for symmetric LFs with jump length distribution 
A(a;) ~ 1x1"^"" (0 < a < 2) the distribution of leapover 
lengths across x — d is distributed like pi{£) ~ £~^~°'/^, 
i.e., it is much broader than the original jump length dis- 
tribution. In contrast, for one-sided LFs the scaling of 
Pi{i) bears the same index a. 

For processes subject to a narrow jump length dis- 
tribution with finite second moment x^X{x)dx the 
crossing of a given threshold value d is identical to the 
first arrival sd x — d 2\. This is no longer true for 
LFs: Intuitively, a particle, whose jump lengths are dis- 
tributed according to the symmetric long-tailed distribu- 
tion X{x) ~ |a;|^^^" (0 < a < 2) is likely to criss-cross 
the point x = d multiple times before it eventually hits 
it, causing the first arrival at d to be slower than its 
first passage across d A measure for the ability to 
criss-cross d is the distribution of leapover lengths, pi{i)- 
Information on the leapover behavior is therefore impor- 
tant to the understanding of how far proteins searching 
for their specific binding site along DNA overshoot their 
target [5], climatic forcing visible in ice core records ex- 
ceeds a given value [6|, or to define better stock market 
strategies determining when to buy or sell a certain share 
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Figure 1: Schematic of the leapover problem: the random 
walker starts at x = and after a number of jumps crosses 
the point x — d, overshooting it by a distance i. 

instead of fixing a threshold price '7| . The quantification 
of leapovers is vital to estimate how far diseases would 
spread once a carrier of that disease crosses a certain 
border Leapover statistics of one-sided LFs provide 
an interesting alternative interpretation of the distribu- 
tion of the first waiting time in ageing continuous time 
random walks just to name a few examples. 

The master equation for a Markovian diffusion process. 
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[X{x - x')P{x\t) 

~\{x' - x)P{x,t)]dx' (1) 

accounts for the influx of probability to position x, and 
the outflux away from x, where A(a;) is a general, normal- 
ized jump length distribution. The time scale for single 
jumps is T. The solution to Eq. ([1]) in Fourier space is 
P{k,t) — e^I^^'^^'^^l*/'^, denoting the Fourier transform 
f{k) = J^^e^^^ f{x)dx by explicit dependence on the 
wave number k. For instance, for the symmetric jump 
length distribution \{x) ~ f7"|a;|~^^", one finds 

P(fc,t) = e-^'°'l^l°* (2) 

with K^"") — cr°' /t, the characteristic function of a sym- 
metric Levy stable law as obtained from continuous time 
random walk theory in the diffusion limit or from the 
equivalent space fractional diffusion equation [lo| . 

In the following we study processes with the long-tailed 
composite jump length distribution 



X{x)/t = e{\x\ - e) [ciQ{-x) + C2Q{x)] /\x 
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where Q{x) is the Heaviside function. For ci = C2, A(a;) 
defines a symmetric LF, and for ci = and C2 > a com- 
pletely asymmetric (one-sided) LF permitting exclusively 
forward jumps. The cutoff s excludes the singularity at 
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X = 0, but can be taken to be small, e ^ 

In the theory of homogeneous random processes with 
independent jumps there exists a theorem, which pro- 
vides an exact expression for the joint PDF p(t, £) of 
first passage time r and leapover length i {£ > 0) across 
X — d for a particle initially seeded at a: = [i3,[ll. We 
here evaluate this theorem, that appears to have been 
widely overlooked, and derive a number of new analytic 
results for p/(t) and pi{t) of symmetric and one-sided 
LFs. With the probability to jump longer than a;. 



M{x) = / \{x')dx\ a: > 0, 



(4) 



the theorem states that the double Laplace transform of 
the joint PDF [H, [11] 
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is given in terms of the multiple integral 

P{u,N = ^-q+{u,d) / ds 

u Jq as 

. ds' Jo 
xM{d + s" - s' - s)ds". 



(6) 



Here, we use the two auxiliary measures q± {u, x) defined 
through Fourier transforms 



q±{u, k) 



= exp < ± 
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and the condition g±(u, 0) = 0. They are related to the 
cumulative distributions of the maximum, Q^{t,d) = 
Pr {maxo<r<t a;(T) < 0?}, and minimum, Q^{t,d) = 
Pr {mino<r<t a;(T) < d}, of the position x{t) such that 
q±{u,d) = u e~'^*Q±{t,d)dt. The complicated inte- 
grals above reduce to elegant results for symmetric and 
one-sided LFs, as we show now. 

For symmetric LFs (ci = C2 = c), the propagator is 
defined by the characteristic function ([2]) with generalized 
diffusion coefhcient if*^"^ = 2cr(l — a) cos(7ra/2)/a. In 
the limit u (long time limit), we obtain from Eq. ([7]) 



q+{u,k) 



,1/2 



VX(^|fc|"/2 
Inverse Fourier transform yields 
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Figure 2: First passage time density p/(r) for symmetric LFs 
with Sparre Andersen universality p/(r) ~ r~"^/^. The curves 
for Q — 0.75 and 1.25 are multiplied by a factor 10 and 100, 
for better visibility. Lines: theory. Symbols: simulations. 



such that from Pf(u) 



- q+{u, d) we find 
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aV7ri^(")r(a/2) 



(10) 



This is the exact asymptotic first passage time PDF of 
symmetric LFs. Fig. [2] demonstrates good agreement 
with simulations results, for which the algorithm from 
Ref. ^lJ| was used to obtain random numbers distributed 
according to Levy stable laws. We note that previously 
only the r"'^/^ scaling was known from simulations and 
application of the Sparre Andersen theorem [^. 
For symmetric LFs, for < a < 2 we obtain that 



M{x) 



K(o') 



2r(l-a) cos(7ra/2) 



x>0. 



(11) 



Using that for symmetric LFs q-{T,x) = g+(r, — x) it 
turns out after some transformations from Eq. ^ that 



_^,sin(^a/2) jd/ir/^ 



TT d + £ 

from which it follows immediately that 
sin(7ra/2) (i"/^ 



TT e°'/^{d + i)'' 



d£, (12) 



(13) 



see Fig.[3l Note that pi is normalized. In the limit a — > 2, 
pi{£) tends to zero if ^ 7^ and to infinity at ^ = 
corresponding to the absence of leapovers in the Gaussian 
continuum limit. However, for < a < 2 the leapover 
PDF follows an asymptotic power-law with index a/2, 
and is thus broader than the original jump length PDF 
A(x) with index a. This is a remarkable finding: while A 
for 1 < a < 2 has a finite characteristic length {\x\), the 
mean leapover length diverges. 
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Figure 3: Leapover density pi{£) for symmetric LPs. 

Consider now one-sided LFs with ci = in Eq. ([3]). In 
this case, the PDF has the characteristic function 

P (k, t) ^ cxp^^- K^°'h\k\°' l-isign(A:)tan(^) |, 

(14) 

where X^"' — C2T{1 — a) cos(7ra/2) /a and A4{x) for x > 
is twice the expression in Eq. pTjl . Eq. ([7]) leads to 

q+{u,k) = ^, C = i^(")Hfc)7cos(^), (15) 

as {—ik)" = [— isign(fc)|A:|]" = exp[— isign(fc)7ra/2]. 
From this we calculate that the Fourier transform of 
(exp(— ur)) = exp(— MT)py (r)c?T can be written as 



i-zky 



{-ik)°' + ucos(7ra/2)/i4:(") 
and we change the variable ik — > — s to find [isj 



.cos(^) d° 



(16) 
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Here, we used the definition of the Mittag-Lefflcr function 



Ea {~ex°')e-''''dx 



(18) 



whose series expansion and asymptotic behavior are [l( 
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^ (-l)"z-i-" 



-;r(l + an) ^ r(l-a[n+l])- 

n— n— 



(19) 



From the relation between and the M^-function [iq . 

e""*M„(t)dt £;a(-u), < a < 1, (20) 
the following result for the first passage time PDF yields 



P/(t) 



cos(a7r/2) d"^^" Vcos (Q!7r/2) 



(21) 



Figure 4: First passage density for one- sided LF {K^"^ = 1). 
The thick lines represent numerical evaluations of the exact 
analytic expression, while the thin dashed lines correspond to 
the asymptotic behavior (|23p . Symbols: simulations. 



The M^-function has the series representation and 
asymptotic behavior with exponential decay 
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1 - a 



a 



(az) 



1/(1-") 



(22) 
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The moments of the il/a-function are obtained through 

z'^MJz)dz = lim(-l)" £^^(-5) = ; 

^ ' ' ds"^ ^ ' r(l + an)' 

(24) 

from which we calculate the mean first passage time 

d"cos(7ra/2) 
" A'(")r(l + a)' 

that is finite and grows with the ath power of the dis- 
tance d. For a — 1/2, we recover the exact form 



(25) 



P^(t) =X("\/ — exp 
Trd 



2d 



(26) 



The first passage PDF Pf{T) is displayed in Fig. 2] in nice 
agreement with the simulations. Note that for a < 1/2 
the tail of X{x) is so long that it is most likely to cross 
X ^ d in the first jump, while for a > 1/2, P/(t) has a 
maximim at finite r > 0. 

To obtain the leapover statistics for the one-sided LF, 
we first note that since P{x < 0,t) = (only forward 
steps are permitted) we have q^{u,k) — 1, and thus 
dq-{u,x)/dx — S{x). Combining Eqs. ([6|) and (O, 

(e-f''') = 1 - lim ^ / / e'l^^'Mid + s' - s) 
^^ouJo Jo 

dq+{u,s 



ds 



ds'ds. (27) 
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Figure 5: Leapover distribution for one-sided LF with d = 10. 



With the small u expansion of the Mittag-Leffler func- 
tion, Eqs. dni) and ^ produce 



dq+{u,x) MC0s(7ra/2) 
Eqs. ([13 and (12S1) inserted into Eq. (gT]) then yield 



sin(7ra) 



-fj.e_ 



leading to the leapover PDF 
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sin(7rQ!) 



(28) 
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which corresponds to the result obtained in Ref. 17| from 
a different method. Thus, for the one-sided LF, the scal- 
ing of the leapover is exactly the same as for the jump 
length distribution, namely, with exponent a. 

The leapover distribution (I30p also provides a new as- 
pect to the first waiting time in a renewal process with 
broad waiting time distribution tp{t) ~ t~^~^ {0 < (3 < 
1) . Interprete the position x as time and the jump lengths 
drawn from the one-sided A(a;) as waiting times t. Con- 
sider an experiment, starting at time to, on a system pre- 
pared at time (corresponding to position a; = 0). Then 
the first recorded waiting time ti of the system will be 
distributed like pi(ii) — tt~^ sm{Tra)tQ /[t"{tn + ti)], as 
obtained from a different reasoning in Ref. [9[. We note 
that the first passage time r in this analogy corresponds 
to the number of waiting events. 

While for symmetric LFs it was previously established 
that the first passage time distribution follows the uni- 
versal Sparre Andersen asymptotics p/(t) ~ t""^/^, here 
we derived the prefactor of this law, in particular, its 
dependence on the generalized diffusion coefficient K^°'\ 



For the same case, we derived the leapover distribution 
Pi{e), that is surprising for two reasons: first, pi{i) is 
independent of K^°'\ synonymous to the noise strength; 
and second, its power-law exponent is a/2, and thus pi{i) 
is broader than the original jump length distribution. 

For one-sided LFs, we recovered the previously re- 
ported leapover distribution and derived the so far 
unknown first passage time distribution. While the 
leapovers follow the same asymptotic scaling pi{£) ~ 
t^^^" as the jump lengths A(x), once more independent 
of K^°'\ the first passage times are narrowly distributed. 
We also drew an analogy between the leapovers and the 
first waiting time in a subdiffusive renewal process. For 
both symmetric and one-sided LFs, extensive simulations 
showed nice agreement with the theoretical results, with- 
out adjustable parameters. 
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